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Abstract — Spectrum sensing is one of the enabling func- 
tionalities for cognitive radio (CR) systems to operate in the 
spectrum white space. To protect the primary incumbent users 
from interference, the CR is required to detect incumbent 
signals at very low signal-to-noise ratio (SNR). In this paper, 
we present a spectrum sensing technique based on correlating 
spectra for detection of television (TV) broadcasting signals. The 
basic strategy is to correlate the periodogram of the received 
signal with the a priori known spectral features of the primary 
signal. We show that according to the Neyman-Pearson criterion, 
this spectra correlation-based sensing technique is asymptotically 
optimal at very low SNR and with a large sensing time. From the 
system design perspective, we analyze the effect of the spectral 
features on the spectrum sensing performance. Through the 
optimization analysis, we obtain useful insights on how to choose 
effective spectral features to achieve reliable sensing. Simulation 
results show that the proposed sensing technique can reliably 
detect analog and digital TV signals at SNR as low as 20 dB. 

I. Introduction 

Due to the increasing proliferation of wireless devices and 
services, the traditional static spectrum allocation policy be- 
comes inefficient. The Federal Communications Commission 
(FCC) has recently opened the TV bands for cognitive radio 
devices, which can continuously sense the spectral environ- 
ment, dynamically identify unused spectral segments, and 
then operate in these white spaces without causing harmful 
interference to the incumbent communication services [1]. 
The IEEE 802.22 Wireless Regional Area Network (WRAN) 
working group is developing a CR-based air interface standard 
for unlicensed operation in the unused TV bands [2]. 

Spectrum sensing to detect the presence of primary signals 
is one of the most important functionalities of CRs. To avoid 
causing harmful interference to the incumbent users, FCC re- 
quires that unlicensed CR devices operating in the unused TV 
bands detect TV and wireless microphone signals at a power 
level of —114 dBm [1]. For a noise floor around —96 dBm in 
the receiver circuitry (with respect to 6 MHz bandwidth and 
a 10 dB noise figure), spectrum sensing algorithms need to 
reliably detect incumbent TV signals at a very low SNR of at 
least —18 dB. This requirement poses new challenges to the 
design of CR systems since traditional detection techniques 
such as energy detection and matched filtering are no longer 
applicable in the very low SNR region [2]. 

In general, there are three signal detection approaches for 
spectrum sensing: energy detection, matched filtering (coher- 



ent detection), and feature detection. If only the local noise 
power is known, the energy detector is optimal [3]. If a 
deterministic pattern (e.g., pilot, preamble, or training se- 
quence) of primary signals is known, then the optimal detector 
usually applies a matched filtering structure to maximize the 
probability of detection. Depending on the available a priori 
information about the primary signal, one may choose one of 
the above approaches for spectrum sensing in CR networks. 
However, energy detection and matched filtering approaches 
are not applicable to detecting weak signals at very low SNR. 
At very low SNR, the energy detector suffers from noise 
uncertainty and the matched filter experiences the problem 
of lost synchronization. To improve sensing reliability, most 
previous studies have focused on the development of coop- 
erative sensing schemes using multiple CRs [4] [5] [6]. An 
alternative approach is to use feature detection provided that 
some information is a priori known. Cyclostationary detection 
exploiting the periodicity in the modulated schemes [7] is 
such an example but requires high computational complexity. 
Recently, Zeng and Liang developed an eigenvalue based 
algorithm using the ratio of the maximum and minimum 
eigenvalues of the sample co variance matrix [8]. 



In this paper, we develop a feature detection-based spectrum 
sensing technique for a single CR to meet the FCC sensing 
requirement. The basic strategy is to correlate the periodogram 
of the received signal with the selected spectral features 
of a particular TV transmission scheme, either the national 
television system committee (NTSC) scheme or the advanced 
television standard committee (ATSC) scheme, and then to 
examine the correlation for decision making. By utilizing the 
asymptotic properties of Toeplitz matrices [9], we show that 
for certain signal models the spectra correlation-based detector 
is asymptotically equivalent to the likelihood ratio test (LRT) 
at very low SNR. In addition, we analyze how the spectral 
features can affect the sensing performance. Specifically, we 
formulate the sensing problem into an optimization problem. 
By solving this problem, we obtain useful insights on how to 
select or design effective spectral features to achieve reliable 
sensing. Extensive simulation results show that the proposed 
sensing technique can reliably detect TV signals from additive 
white Gaussian noise (AWGN) at SNR as low as -20 dB. 
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(a) The measured NTSC channel spectrum in UHF 
Channel 51 (San Diego, CA, USA). 



ATSC Signal (CH 19) 




(b) The measured ATSC channel spectrum in UHF 
Channel 19 (San Diego, CA, USA). 

Fig. 1. The estimated power spectra in NTSC and ATSC channels. 



II. Spectra Correlation Based Spectrum Sensing 

Before presenting the spectrum sensing technique, we first 
briefly review the TV transmission schemes. 

A. TV Signal Characteristics 

A typical TV channel occupies a total bandwidth of 6 MHz 
and its power spectrum density (PSD) describes how the signal 
power is distributed in the frequency domain. Fig. [T] (a) and 
(b) illustrate the PSD functions of both NTSC and ATSC 
signals. NTSC is the standardized analog video system used 
in North America and most of South America. The power 
spectrum of an NTSC signal consists of three peaks across 
the 6 MHz channel, which correspond to the video, color, 
and audio carriers, respectively. On the other hand, ATSC is 
designed for the digital television (DTV) transmission, and 
it delivers a Moving Picture Experts Group (MPEG)-2 video 
stream of up to 19.39 Mbps. The ATSC spectrum is relatively 
flat but has a pilot located in 310 kHz above the lower edge 
of the channel. 

We find that both NTSC and ATSC signals have distinct 
spectral features, which are constant during the transmissions. 
This observation motivates us to design a spectrum sensing 
technique for TV signals by exploiting these a priori known 
spectral features. 
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B. Sensing Strategy 

The spectrum sensing problem can be modeled into a binary 
hypothesis test at the Z-th time instant as follows: 

H a : y(l) = v(l), 1 = 0,1,2,...; 
Hi : y(l) = x[l)+v(l), 1 = 0,1,2, 

where y(l) is the received signal by a secondary user, x(l) 
denotes the transmitted incumbent signal, and v(l) is as- 
sumed to be complex zero-mean additive white Gaussian 
noise (AWGN), i.e., v(l) ~ £/V(0, cr^). We assume that the 
signal and noise are independent. Accordingly, the PSD of the 
received signal <Sy(w) for different hypotheses can be written 
as 
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Hi : S y (uj) = S x (u) + <, < u < 2tt, 



(2) 



where Sx(u) is the PSD function of the transmitted primary 
signal. Our objective is to distinguish between Ho and Hi by 
exploiting the unique spectral signature exhibited in 5*x(w). 

Generally, we can obtain an estimate of the PSD of the 
observations through various spectral estimation algorithms, 
and here we focus on the periodogram, i.e., the squared 
magnitudes of the n-point discrete-time Fourier transform 
(DFT) of the ri-point received signal, denoted 



S { y\k), fc = 0,l, 
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On the other hand, we suppose that the n-point sampled 
PSD of the signal under detection, S^\k) = S x {2-Kk/n), 
is known a priori at the receiver. To detect the presence of a 
TV (NTSC or ATSC) signal, we perform the following test: 
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where 7 is the decision threshold. Namely, if the spectra 
correlation between S^\k) and Sy (k) is greater than the 
threshold then we would decide Hi, i.e., presence of the signal 
of interest; otherwise, we would decide Ho, i.e., absence of 
the primary signal. 

III. Asymptotic Optimality 

In this section, we show that the proposed spectrum sensing 
technique © is asymptotically optimal at very low SNR in the 
Neyman-Pearson sense. The asymptotic optimality is in the 
sense that, as shown in Theorem 1 below, the decision statistic 
T n asymptotically approaches the likelihood ratio decision 
statistic for low SNR and large observation length. 

A. LRT at Very Low SNR 

Considering a sensing interval of n samples, we can rep- 
resent the received signal and the primary transmitted signal 
in vector form as y = [y{0), y(l), ■ ■ ■ ,y(n — 1)] T and x = 
[x(0), x(l), . . . , x(n — 1)] T . Since TV signals are perturbed 
by propagation along multiple paths, it may be reasonable to 
approximately model them as being a second-order stationary 
zero-mean Gaussian stochastic process, i.e, 

x^CAA(0,S„) (5) 



where 
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is the covariance matrix. Consequently, (Q~|i is equivalent to 
the following hypothesis testing problem in the n-dimensional 
complex space C n : 



H : y~ OA/" (0 jO ^l) 
HH y~CjV(0,S n + (72l) 



(7) 



where I is the identity matrix. The logarithm of the likelihood 
ratio is given by [10]: 



logL(y) =2nloga v - logdet (£„ + o%l) 
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Incorporating the constant terms into the threshold, we obtain 
the logarithmic LRT detector in the quadratic form as 
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which is the optimal detection scheme according to the 
Neyman-Pearson criterion. This detector is also known as a 
quadratic detector. 

From the Taylor series expansion, we have 
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where the convergence of the series is obtained if the eigen- 
values of (T ! 7 2 S n are less than unity. This condition always 
holds in the low SNR regime where al grows sufficiently 
large. For weak signal detection in the very low SNR region, 
i.e., det 1 /" (S„) <§; ct, 2 ,, ( fT0T > can be approximated as 



(a 2 J- 
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Plugging ( fTTb into (O, we obtain 
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Hence, the optimal LRT detector at very low SNR is given by 
i H\ 



Tlrt,-, 
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where 7lrt = '/n. 

B. Asymptotic Equivalence 

Now we show that our proposed spectra correlation-based 
detector (0]i is asymptotically equivalent to the LRT detector 
at very low SNR ( TT2l . Consider a sequence of optimal LRT 
detectors as defined in ( flZb 
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7lrt,?j = ^-y T S n y ^ 7lrt, n = 1,2, 
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Likewise, we define a sequence of spectra correlation detectors 



1 n— 1 

r«=-y)S , i? ) (fc)4 T,) (*). n = l,2,.... (14) 

fc=0 

Note that the LRT detectors are performed in time domain 
while the spectra correlation detectors are in frequency do- 
main. The asymptotic equivalence of these two sequences of 
detectors is established in the following theorem. 

Theorem 1: The sequence of spectra correlation detectors 
{T n } defined in (fl4l are asymptotically equivalent to the 
sequence of optimal LRT detectors {Tlrt,™} at very low SNR 
defined in ( fl3l ). i.e., 

lim |T LRT ,„ - T n \ = 0. (15) 

n — >oo 

Proof: The proof is sketched in Appendix |A] ■ 

IV. Spectral Feature Selection 

In this section, we study the effect of spectral features on 
the detection performance. Although the sensing algorithm 
cannot control or change the spectral features of transmitted 
signals since these features are completely determined by 
the incumbent transmitter, we can obtain through analysis 
important insights to identify the best features for the signal 
detection. These insights are also useful for system engineers 
to design robust signals that can be reliably detected at very 
low SNR. 

We first consider the case where there is no primary signal 
in the band of interest. Under hypothesis Ho, we have 
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is the average power transmitted across the whole bandwidth. 
On the other hand, by exploiting the fact that the periodogram 
is an asymptotically unbiased estimate of the PSD [11], we 
have for sufficiently large n, 

-. n — 1 

E[T„, 1 ] = -^E[^" ) (fc)]5^ ) (fc) 
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1 n—1 



(18) 



under hypothesis Hi. 

Here we shall use the difference between E[T„ o] and 
E[T„.i] to determine the detection performance. Suppose that 
we can control the spectral mask {S^\k)} of the transmitted 
signal, we would like to maximize the difference between 
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Fig. 2. A geometric illustration of the non-convex optimization problem 
formulated in 1201 . 

E[T n , ] and E[T nA ], i.e., 

maximize E[T ni i] — E[T„.o] 

s.t. ±J2sP(k) = P x 

S ( x\k)>0, fc = 0,l,...,n-l (19) 

with the optimization variables {S^ (fe)}fc=o' For large n, 
this problem is equivalent to 

ri — 1 2 

maximize [^x ' (^) 

fc=0 
n-l 

s.t. ^5i n) (fc) = nP x 

s£°(&)>0, fc = 0,l,...,n-l (20) 

which maximizes a convex function over a hyperplane. In the 
sequel, we will show how to solve this nonconvex optimization 
problem. 

To solve (l20l . we first look at its geometrical representation, 
as shown in Fig. |2] It is easy to see that optimal solutions fall 
into the intersection of the convex surface of the objective 
function and the hyperplane determined by the constraints. 
Thus, the optimal solutions are given by 

f S { p(j) = nP x , je {0,1,..., n-l} 

[ 5^ ) (fc) = 0, < k < n - 1, and k^j 

for any arbitrary j, implying that all the transmit power 
is concentrated in a single frequency bin. Accordingly, the 
optimal value of (l20l i is 

E[T„,i] - E[T n , ] = n 2 P*. (22) 

On the other hand, the worst case occurs when 

S$\k)=P a , fe = Q,l,...,n-l (23) 

which makes E[T ni i] — K[T n fi] = nP 2 . In this extreme case, 
the spectral mask function is flat across the spectrum. A repre- 
sentative example of such a spectral feature is the white noise- 
like signal. This result is consistent with our intuition since it is 



TABLE I 

ITU Pedestrian-B Multipath Channel - Power Delay Profile 
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Fig. 3. The missed detection rate of the proposed spectrum sensing algorithm 
for ATSC signals, with a false alarm rate less than 0.001. The detection interval 
is 24 ms. 

generally difficult to distinguish a white Gaussian signal from 
additive white Gaussian noise. The optimal detector for such 
a case is the energy detector (radiometer) [10] [3] provided 
that the noise power is perfectly known. 

V. Simulation Results 

In this section, we numerically evaluate the proposed spectra 
feature correlation-based spectrum sensing algorithm for both 
ATSC and NTSC signals. 

First, we obtain the "clean" baseband TV signals by captur- 
ing the TV signals in the RF front-end and then transforming 
the signals from the ultra-high frequency (UHF) bands to the 
baseband. Through processing, the TV signals in the baseband 
are sampled at a rate 6 x 10 6 samples/sec, with 6 MHz 
bandwidth. The TV data are divided into a number of blocks, 
each of which is 6 msecs. The spectral features are obtained 
by computing the (averaged) periodograms of the clean TV 
signals. 

We study the sensing performance for two channel models: 
AWGN and multipath fading. For the multipath fading channel 
model, we apply the ITU Pedestrian B model, whose power 
delay profile is given in Table U The root mean square (RMS) 
delay spread of the ITU Ped-B model is 633 ns. We pass the 
clean TV signals through the channel models to simulate the 
real wireless propagation environment. 

For both ATSC and NTSC signals, we choose the test 
thresholds such that their false alarm rates are less than 0.001. 
We use the white Gaussian noise to test the sensing algorithms 
to make sure the false alarm rates are less than 0.001. Once we 
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Fig. 4. The missed detection rate of the proposed spectrum sensing algorithm 
for NTSC signals, with a false alarm rate less than 0.001. The results for 
AWGN channels are not shown because they are lower than 10 — 4 in the 
above SNR regions. The detection interval is 24 ms. 



find the test thresholds for ATSC and NTSC signals, we can 
simulate and calculate the missed detection rates. For each 
SNR value, we simulate the sensing algorithms for 25,000 
realizations. The simulation results for ATSC and NTSC are 
plotted in Figs. [3] and |4] 

For both ATSC and NTSC signals, the spectral feature 
detector can reliably detect the signals at SNR= —20 dB with 
a missed detection rate less than 0.01 in the AWGN channels. 
It can be observed that the NTSC signal is easier to detect 
than the ATSC signal since the NTSC signal has three sharp 
spectral features while the ATSC contains a large amount 
of flat spectrum and has only one feature corresponding to 
the pilot. This observation is consistent with our optimization 
analysis in Section ITVl 



VI. Conclusion 



In this paper, we have proposed a spectral feature detector 
for spectrum sensing in CR networks. The basic strategy is to 
use the correlation between the periodogram of the received 
signal and the a priori spectral features. Using the asymptotic 
properties of Toeplitz and circular matrices, we have shown 
that this spectral feature detector is asymptotically optimal at 
very low SNR and with a large block size. In addition, we 
have performed optimization analysis on the effects of spectral 
features on the sensing performance. The analytical results 
show that the signals with sharp spectral features are easier 
to detect compared with those with relatively flat spectra. 
The simulation results show that the proposed spectral feature 
correlation detector can reliably detect analog and digital TV 
signals at SNR as low as -20 dB. 



Appendix A 
Proof of TheoremQ] 

From the definition of the two test statistics, we have 

Hm |T LRT ,„-T„|= lim - |y*£„y - y*W*\ n W n y\ 

n— >oo n — >oo fi 



where 
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is a diagonal matrix with the PSD of the incumbent signal in 
the diagonal, and W n is the DFT matrix defined as 



W n = 



- 1 


1 


l 


1 


Wn 


2 
11)- 


1 


o 

W 

a 


4 

W 

n 



., n — 1 2(n—l) 

1 1 Wn- ' 



1 

w n-l 
2(n-l) 



(n-l)(n-l) 
Wn 



(26) 



with w n = e 3 2n / n being a primitive nth root of unity. 
Consequently, 

lim |T LKr , n - T n | = lim - |y* (S n - W*A n W n ) y| 

n — >oo n — *oo fl 

= lim - |y* (£„ - C n ) y| (27) 



where C„ = W*A n W n is a circular matrix. It has been 
shown in [9] that the Toepliz matrix S n is asymptotically 
equivalent to the circular matrix C n since the weak norm 
(Hilbert-Schmidt norm) of S„ — C n goes to zero [9], i.e., 



lim ||S n -C„|| =0. 

n — >oo 

Thus, we can establish (l5i from ( 1271 ). 
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